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Abstract. We study the energy of relativistic electrons and positrons inter- 
acting via the second quantized Coulomb potential in the field of a nucleus of 
charge Z within the Hartree-Fock approximation. We show that the associated 
functional has a minimizer. In addition, all minimizers are purely electronic 
states, they are projections, and fulfill the no-pair Dirac-Fock equations. 



1. Introduction 

Heavy atoms should be described by relativistic quantum mechanics. It is com- 
monly believed quantum electrodynamics (QED) yields such a description. For- 
mally the Hamiltonian is given (Bjorken and Drell Formula 15.28]) as 

H = Jdx: ip*(x) \a ■ (iV - y/aA(x)) + mf3 - ^1 ij)(x) : 



(1) 

+ 2J""J" S |x-y 



Jdx J dy^M^rim^ + _L J r3 . B (x) 2 + A (x) 2 : dx 



where the normal ordering denoted by colons is with respect to a given choice of 
the one-electron space. However, it is not clear how this expression can be self- 
adjointly realized as a positive operator. To simplify matters we omit the energy 
of the transverse radiation field coupled to the current, i.e., we set B = and 
A = in the above expression. We will also regularize the Coulomb interaction 
of the electron-positron field by normal ordering it completely. Both assumptions 
are simplifications. The first one can be justified by the physical wisdom that the 
presence of self-generated magnetic field is physically known to be small compared 
to the relativistic effects in heavy atoms. Moreover, the inclusion poses serious 
technical problems that we presently cannot solve. The second assumption ignores 
the vacuum polarization effects which are also small comared to the relativistic 
effects and would contribute to the Lamb shift only. 

Based on an interesting observation of Chaix, Iracane, and Lions |HlE3i Bach 
et al. Ej showed positivity for the corresponding quadratic form without any 
constraints on the charge of the state, if the one-electron subspace is appropriately 
(Furry picture) chosen. In particular they showed that the vacuum state (particle 
number equal to zero) has energy zero. 
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In order to do so, they proved that the positivity of H on generalized Hartree- 
Fock states (quasifree states with finite free kinetic energy) is equivalent to the pos- 
itivity of the Hartree-Fock functional (see the Definition given in H2.IJI). a functional 
on density matrices 7, where the charge is given as tr 7, which thus is unrestricted. 

For describing atoms one needs, however, to restrict to states with prescribed 
charge q. In order to implement this problem, we subtract from the energy the rest 
mass m, which - as we will see - will allow us to relax the constraint tr 7 = q to 
tr 7 < q. 

Given an electron space by the positive spectral subspace of a Dirac operator 
with a mean-field potential, we show that there exists a minimizer of the associated 
Hartree-Fock functional in a suitable set, if q < Z (Theorem I3.9|) . Moreover, if q 
is a positive integer, all minimizers are projections of maximal rank. In particu- 
lar, minimizers are purely electronic and are projections (Theorem I4.3j) . Finally, 
we show that the eigenfunctions of any minimizer 7 fulfill the no-pair Dirac-Fock 
equations with the q first (positive) eigenvalues (Theorem 14. 5fl . 

According to Mittleman ^7], the most stable, i.e., highest ground state energy 
is the physical ground state. Thus one should maximize among the allowed one- 
electron subspaces, yielding a max-min variational problem. The resulting Euler 
equation should, on a heuristic level, give the Dirac-Fock equations which were 
treated by Esteban, Sere, and Paturel [I'M 1191 IT2] . There are indications that this 
latter question might be answered affirmatively only under additional hypotheses: 
In a recent work [o] an atom with total charge q— Z < is considered. If the ground 
state of the noninteracting problem with N — q electrons corresponds to a closed 
shell, then for small interaction maximizing over one-electron subspace yields the 
Dirac-Fock equations in the non-relativistic limit. However, if the noninteracting 
problem corresponds to an open shell, then, in the same limit, the max-min pro- 
cedure does not yield a solution of the Dirac-Fock equations with self-consistent 
projector as considered by Esteban and Sere [T^J. While this result is perturbative, 
it indicates on the one hand side that the Dirac-Fock equations and the Mittleman 
principle might agree in the case of filled shells whereas in the unfilled shell case it 
might give different results which raises the question which procedure is physically 
relevant, a problem that we have to leave open at this point. 

We add a short guide through the paper for the orientation of the reader: Sec- 
tion contains some basic material. We define the set of density matrices that will 
be allowed. There will be two types of density matrices, the charge density matrices 
7 for the electron-positron field and the density matrices S giving the screening of 
the one-particle Dirac operator that defines the electron subspace. In addition this 
section contains some basic estimates on the direct and exchange energy. Section [21 
contains the actual minimization. We first show that the elimination of positrons 
lowers the energy (Lemma 18.1(1 : next we show that the density matrices can be 
restricted to finite rank (Lemma I3.3J1 . and the minimization under the constraint 
tr 7 < q < Z gives a minimizer - if existing - with charge equal to q (Lemma 13. 711 . 
This allows us to show the existence of minimizers (Theorem 13. 9fl . In Section 0] we 
investigate the minimizers. They turn out to be projections that fulfill the no-pair 
Dirac-Fock equations (Theorem I4.5JI . In the last section we give an outlook with 
respect to the above mentioned program of Mittleman. We derive the correspond- 
ing Euler equation (Theorem 15. 1(1 . However, we are not able to show that there is 
a maximizer. 
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2. Definition of the Problem 

A single relativistic electron or positron in the field of a nucleus of charge Z can 
be described by the Coulomb-Dirac operator 

1 Z 

(2) D z ■= a. ■ -V + 771/3 - a-r—r , 

i |x| 

where a is the Sommerfeld fine structure constant. 

The operator Dz is self-adjointly realized in Sj := L 2 (R 3 ) (gi C 4 and essentially 
self-adjoint on C£°(M 3 \ {0}) ® C 4 , if aZ G (-\/3/2, \/3/2). Here, we will even 
assume 

(3) aZ £ [0, V3/2) and a > , 
throughout the paper. The domain of Dz is H 1 (G) where 

G :=R 3 x {1,2,3,4} . 

(Landgren and Rejto ^5]> Theorem 2.1). For Z = 0, Z?o is just the free Dirac 
operator. 

In the present paper, we consider a larger class of Dirac operators, namely Dirac- 
Fock operators. They are Hamiltonians for a relativistic particle in a mean field 
created by other particles. For that purpose, we will consider operators with an 
additional mean field potential 

(4) wW=(pW-XW, 

where tp^ and X^ will be defined below. 

Forp G [l,oo), we denote by 6 p (Sj) B(Sj) | tr \A\p < oo}, and by Soo(£) 

the space of compact operators on $). 

Definition 2.1. We denote by F the Banach space of all self-adjoint operators 8 
on S) such that the norm \\S\\f '■= tr ( l-Dol 1 ^ 2 \$\ l-Dol 1 ^ 2 ) is finite. 

For a given element 8 £ F, we denote by (A„) the sequence of its eigenvalues and 
by (£„) a corresponding orthonormal basis of eigenvectors; the associated integral 
kernel d(x, y) is 

(5) S( X >y) '■= X! X n£,n(x)^n(y) ■ 

n 

(It is convenient to introduce the notation x — (x, s) for an element of G and dx 
for the product of the Lebesgue measure dx on R 3 with the counting measure in 
{1, 2, 3, 4}.) Associated with 8 is its one-particle density 

4 

(6) P5 (x):=^^A„|C„(x)| 2 , 

s— 1 n 

its electric potential 

(7) ^(x)=/ ^Uy, 

Jn* |x - y| 

and its exchange operator X^ s ' 

(8) -^-dy . 

The difference of these two operators is the mean field potential W^ s ' defined in 
(0J. Next, we define for the given 8 the Coulomb-Dirac operator associated to 8 as 

(9) L> (<5) := D z + aW (s] . 
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As shown in the Appendix 1X1 (Lemmata TA.7I and IA.8|) . the operator is 
bounded implying that flW is self-adjoint with the same domain as the Coulomb- 
Dirac operator Dz which for aZ € [0, v3/2) is identical with the domain of Do- 
Moreover (see Lemma lA.QJI . is relatively compact with respect to Do which 
implies 

(T eS s(-D (6) ) = cr css (D a ) = (-oo, -m] U [m,oo). 

Finally, since l-D^I > c a ,z,m,s\Do\ by definition of c at z,m,s (Equation ||H3) and 
since |£>o| > 0, the operator has a bounded inverse, as soon as c at z, m ,6 > 0. 

The one-electron states are vectors in f)+ = A+Sj where A + is an orthogonal 
projection on Sj, whereas one-positron states are charge conjugated states in A-fi, 
where A_ := 1 — A + ( 22 ). We will take A + to be the projection A^ onto the 
positive spectral subspace of the Dirac-Fock operator D^ s \ i.e., 

(10) A+:=Af :=X[o,oo)P (5) ), 

where xi denotes the characteristic function of the set /. Thus, the choice of 5 fixes 
the definition of the spaces of electrons and positrons. 
The Coulomb scalar product is 



(ID D[ P ,a]:=U dx/ dy-^ ff(y) 



j(x,y)-y'(x,y) 

d v — i i — 



2 7r3 Jrs |x — y 
The exchange scalar product for 7, 7' £ F is 

(12) £[ 7 , 7 ']:=i|dz| 

Lemma 2.2. Assume that 7,7' S F. Then 

(13) WAt-M < J|Mlitr(V=A|y|) , 

(14) |S[7,7']| < D[p| 7 |,P|y|] • 

Proof. Expanding 7 and 7' in their respective bases of eigenfunctions (see ©), wc 
get by the Cauchy-Schwarz inequality 

(15) 



i(x,yMx,y) dxdy 



x-y 



< / / g^M^^™ d ^ = £ ^w^M dxdy , 

which shows that it suffices to prove (I13|) . 

To this end we remark that by Kato's inequality 

(i6) J dz\u*)\ 2 J d y l j^j\ ^ f (C MO • 

The claimed inequality follows now by multiplication with |A^A^| and summation 
over fi and v. □ 

We will also need the following result of Bach et al (0, Equation (30)) : 
Lemma 2.3. IfjEF, then 

£[7,7] < \ tr ( 7 *x/^A 7 
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Given an operator A on the Hilbcrt space Sj the symbols A ++ = AS AAS , 

A+- = A^AA^, A_+ = A^AA^, and = A^AA^ denote the matrix 

elements of the decomposition of A with respect to the splitting of the Hilbert 
space given by aS^ and A_ (We assume that the operators A±* leave the domain 
of A invariant). 

Rclativistic electrons and positrons are described by one-particle charge density 
matrices 7 with certain additional properties reflecting the charge of the particle, 
and the fact that they are Fcrmions and thus obey the Pauli principle. 

Definition 2.4. Given 5 G F and q G M+ we define the following sets of one- 
particle diagonal charge density matrices: 

(17) := {7 G F I - A W < 7 < Af , A<? 7 Af = 0} , 

(18) := {7 G I < tr 7 < q} , 

(19) Sg> := {7 G S<*> I tr 7 = «} . 

We note that all sets are closed subsets of F. Furthermore, the first two are 

convex. Note also that for 7 G S^ s ' we have 7++ > and 7 < which follows 

directly from the definition. We also observe that : 

(20) 7 ++ < 7++, 

(21) 7-- < -7—, 

which permits to get (in the case of diagonal density matrices) 

(22) tr(|L>W|| 7 |) = tr(L>^ 7 ) > tr(|D^| 7 2 ) . 

The elements of are the one-particle (renormalized) charge density matrices 7 
of the electron-positron field. The trace is its total charge. Since we are interested 
in describing atoms we want to fix the charge of the electron-positron field to be q 
and minimize the energy over the set SgJ . For technical reasons we will also use 




We wish to point out that he derivation of the variational spaces of one-particle 
charge density matrices as done in [F] does not give the extra condition A_ 7AS? = 
that appears in the definition of . A formal calculation shows that if we do 
not assume that the one-particle density matrices have off diagonal terms equal 
zero, then in most cases, Lemma 13.11 and Theorem 14.31 and 14.51 do not hold. In 
particular, minimizers will contain electron-positron pairs. 

The projections A^ can be physically interpreted as the one-particle density 

matrix of the Dirac sea. In particular, A_ is the Dirac sea under the influence of 
a nucleus of charge eZ and an electron-positron distribution given by 8. 
For later purposes we also introduce (unrenormalized) density matrices as 

r:=A^+ 7 , 

representing all electrons including those of the Dirac sea. 

The unrenormalized density matrices are nonnegative expressing the fact that 
positrons occur in this picture only as 'holes' in the Dirac sea. 

The energy of a system of electrons and positrons in Hartree-Fock approximation 
is given by the functional 



(23) 



7 h-> tr (Dzj) + aQ[j, 7] , 
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where 

(24) Q ■ SltxSV - R 

1 ' W " (7,7') - Dip^p^-Eiw']. 

As explained above we are primarily interested in the infimum of £\„(.s)', for 

technical reasons we will also consider, for [i £ R, 

-> R 
7 i— » £ (7) — /itr 7 



(25) £ M 



Lemma 2.5. For any /x € R, the energy functional £^ is well defined and contin- 
uous in the \\ ■ \\p norm. 

Proof. The lemma is an immediate consequence of the definition of the norm, 
Lemma f2.3l together with the fact that j 2 < \ j\ and Inequalities 1)13(1 and 1)14(1. □ 

3. Minimization of the Energy 
3.1. Reduction to Electrons. 

Lemma 3.1. Assume q > 0, 7 £ Sg q ; , < S G F. Moreover, assume 

c a ,z,m,6 > 7ra(l/4 + max{tr 5,q}) . 

Then there exists a nonnegative ~f e £ Sg^ and R G Sg^ such that 7 = ~f e + R and 
£{le) < £(l)- I n addition, equality can only occur if < 7. 

Physically speaking this lemma says that it is favorable to have no positrons in 
the system and to restrict the minimization to electron states. 

Proof. Using the spectral decomposition of 7, one can easily construct j e and R 
such that 7 e € Sg^ , R <E SgJ and 7 = 7 e + R. We also note that we can pick R ^= 0, 
if 7 ^ 0, and that we can pick i? = 0, if 7 > 0. We have 

£(7 e ) - 5(7) = - tr (D Zj R) - aQ[i?, i?] - 2aKQ[i?, 7e ] 

(26) = - tr {D {S) R) - aQ [R, R] - 2a$tQ [R, 7e ] + 2a^tQ [R, 5} 

< - tr (D^R) + aE[R, R] - 2dRQ[R, 7J + 2a5RQ[i?, 5] , 

where in the last inequality, we used positivity of D [pn , pa] . 

Let i?+ and i?_ be respectively the positive and negative part of R, i.e. 

R = R+- R- with R± > . 

Using Lemma IA.8I and the positivity of 7 e yields 

KQ[i?+,7 e ] > and thus - UQ[R^ e ] < ^Q[R-, le ] . 

Similarly, we have $IQ[R, S] < 5}. Therefore, from Inequality ((26(1 we obtain 

(27) £(j e )-£(j) < - tr (D 1 -^ R) + aE[R, R] + 2aMQ[R-,"f e ] + 2a%tQ[R + ,S] . 
From Lemma 12.31 and the definition of c a> z,m,S i n <|87(1 . we get 

(28) E[R, R] < - tr (\D \R 2 ) < tr (|X)W \R 2 ) < tr (D^R) , 

where the last inequality is a consequence of l)22[l. 

Using 1(140 ■ I|13|). and, finally, the definition of c Q ,z,m.,<5 in l|S7|) we have 

KQ[fl_,7e] = L>[p ii _ ) p 7e ]-SR J B[ii_ ) 7 e ]<2D[ /9R _,p 7e ] 

(29) <ftr 7 etr(| P |^)<-^^tr(|^)|i?_). 
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With exactly the same arguments as above, exchanging i?_ by R + and j e by S 
we get 

(30) XQ[R + ,S\ < ntlS tr(\D^\R + ). 
Since tr (\D^\\R\) = tr (D^R), we have, using ^ZJ-^JJ 

(31) £ (7e ) - £(7) < f-1 + + tr (U«iQ < , 

under the hypothesis of the theorem. In addition we note that the last inequality 
is strict unless R = or the prefactor vanishes. □ 

3.2. Reduction to Density Matrices with Finite Spectrum. 

Lemma 3.2. Assume < q, < 5 E F and < 7 G Sg. . Then there exists a 

sequence of finite rank density matrices 7^ G Sg 1 such that ~ i\\f — > as 
K — > 00. 

Proof. Let (£,k)kefi be a complete set of eigenfunctions in H 1 (G) of 7 associated 
with the eigenvalues A&. Since 7 > 0, we have 7 = 7A+ and, for all k, we have 
& G A^ 53. 

Assume that 7 is not already of finite rank. Then, since 7 is trace class, there 
exist infinitely many eigenvalues of 7 in (0, 1). Let us pick A„ G (0, 1), one of these 
eigenvalues. 

Set €k '■= q — X)fc=i ^k- Then ek is a nonnegative monotone decreasing sequence 
tending to zero. Define -y K := J2k=i ^k\£k)(£k\ + £id£na)(£no !• 

We assume that K is chosen sufficiently large so that K > no and A„ + ek < 1 • 
Obviously each jk is nonnegative, belongs to and has finite rank. 

We now show that jk — > 7 in F-norm as K — > 00. We have 

00 

1-1K= Afc|ffc}(6i| -f*-|Cno)(£n a | • 

fc=RT+l 

Thus, 

00 

(32) || 7 -7aHIf< ^ A fc tr(|i?o||efc)<&|)+ejf«no,|i5o|U)- 

fc=A+l 

The first in the right hand side tends to zero, since |L>o||7| G 61(53), and the second 
tends to zero, since ek — > 0. □ 

The following is an immediate consequence of the continuity of £ in the F-norm 
and the preceding density result. 

Lemma 3.3. Assume that q > 0. Then 

inf{^( 7 ) I < 7 G 4?} = mf{^(7) I < 7 G 5g, rank( 7 ) < 00}. 

3.3. Reduction to Projection. Following an argument of Bach 0, we get 

Lemma 3.4. Assume 5 e N, 5 E F, < 7 6 Sg^ with finite rank. Then there 
exists a projection P E Sg^ such that 

(33) £{P) < £(7) ■ 
Equality in (|33|l holds only if 7 is already a projection. 
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Proof. Suppose that 7 is not a projection. Then 7 has an eigenvalue A G (0,1); 
we denote a corresponding normalized eigenvector by u and observe that it is in 
H 1 {G). Since tr (7) G N, there exists at least a second eigenvalue fi G (0,1); we 
denote a corresponding normalized eigenvector by v. We set 

(34) 7 := 7 + eS , 
where S 1 := — We get 

(35) £ (7) - £(7) = e(tr (L> Z S) + 23?Q[ 7 , 5]) + e 2 Q[S, 5] . 

By explicit computation and use of the Cauchy-Schwarz inequality, one can show 
that : 

(36) Q[S, S] < , if S is a difference of two orthogonal rank one projections. 

Now - depending on the sign of the coefficient linear in e - we lower the energy by 
increasing or decreasing e from zero, until one of the constraints < A + e,/i — e < 1 
forbids any further increase or decrease of e. This process leaves all the eigenvalues 
of 7 unchanged except for fi and A, one of which becomes either or 1. 

Since there are only finitely many eigenvalues of 7 strictly between zero and 
one (even if they are counted according to their multiplicity) , iterating this process 
eliminates all eigenvalues that are strictly between and 1, i.e., we have found a 
density matrix P such that P 2 = P. □ 

Remark 3.5. Following the same method in the case of q > 0, not necessarily 
integer, we can show that, given < 7 G Sg q , there exists P equals a projection 

plus a rank one operator such that tr P = tr 7, P > and £{P) < £ (7), with 
equality only if 7 is already a projection plus a rank one operator. 

3.4. Lower Bound on the Energy. In this subsection, we show that for suffi- 
ciently small fine structure constant a and atomic number Z ', the energy is bounded 
from below. 

Theorem 3.6. Assume < 8 G F, and c a ,z,m.s > 7ra(l/4 + max{tr 6, q}) > 0. 
Then, for all 7 G S^ S) , £ (7) > 0. 

Proof. By Lemma 13.11 we need to consider only positive 7's. In this case, Ijl4(l 
implies Q[7,7] > 0. Now, for / 6 A^fj, using the definition ij%7|l of c a ^z,m,s, 
Inequality l|13(l and the positivity of E^ s \ we obtain 

(f,D z f) = (f, \D&\f) - a(f, W&f) > c a , z , m , s (f, If I/) ~ ~<*trS(f, \p\f) > . 

Thus, under our hypotheses, £(7) > tr (Dzj) > 0. □ 

3.5. Reduction to Density Matrices of Maximal Charge. 

Lemma 3.7. Assume 0<q,0<6(£F and suppose that, for all < p < q and all 

< 7 G Sp S \ the operator has infinitely many eigenvalues in (0,m). 

Then 

inf {f ro ( 7 )|7 e 4?.7> 0} = inf {f m ( 7 )|7 € ^,7 > o} . 

J/ m addition < 7 is a minimizer of £ rn in Sq S \ it follows that tr 7 = q. 

Proof. That the left side bounds the right side from above is obvious. To prove 
the converse inequality, we assume < 7 G S { q S) , with tr 7 < q. By Lemma EHl 
and Remark 13.51 we can assume that 7 is a projection plus a rank one operator; in 
particular its range is finite dimensional. Since by assumption the discrete spectral 
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subspace X of A^D^A^ is infinite dimensional we can find u G Sj+ nlfl 
with < \\u\\ < 1 and define 7 := 7 + A with A := \u >< u\. We then get 

£ m ( 7 + A) - £ m ( 7 ) = tr (D Z A) + 23tQ( 7 , A) - m\\u\\ 2 = (tt, (£ (7) - m)u) < . 

Thus, this construction yields a density matrix 7 with strictly smaller energy and 
a trace that can be made bigger by min{l, q — tr 7}. Iteration of the construction 
yields the desired result. This proves both claims. □ 

3.6. Existence of a Minimizer. We wish to show the existence of a minimizer 
by weak lower semi-continuity of the functional on a minimizing sequence and weak 
compactness. However, we are faced with the problem that we are minimizing on 
charge density matrices and the fact that the Coulomb potential is not relatively 
compact with respect to the relativistic energy \p\. The first problem has been 
addressed by Solovej [23 m the non-relativistic context. To handle the second 
problem we will decompose the one-particle part of the energy. 

Theorem 3.8. Assume < c a z,m,6> < 5 €z F and q G N. Furthermore, suppose 
that £ m is bounded below on Sg^ and 

inf{£ m ( 7 ) I 7 e S^} = inf{£ m ( 7 ) | 7 G s£\ 7 2 = 7} • 
Then the energy functional E m \ „w has a minimizer. 

Proof. Let j n be a minimizing sequence of orthogonal projections in Sg^ . 

Step 1. Weak Limit of the Minimizing Sequence: Since E m {ln) and <p 5 ' 
are bounded (Lemma IA.7J1 (and thus also X^ s \ by Lemma f A. 81 and W^) there 
exists a constant C such that, for any p > 1 

C > S m (7n) +aq\\W^\\ > tT(\D^\j n ) > c a , z>m , 5 tr (|D |7») = c ,z 1 m,i||7n||f 

> c Q , z , m ,iiii^or /2 7„i^or /2 ii P . 

Thus, according to Banach and Alaoglu, there exists, for all p > 1, 7oo with 
II |-Do| 1 ^ 2 7c3o|-Do| 1 ^ 2 ||p < 00 and a subsequence of 7n - also denoted by 7n - such 
that for all B with \\B\\ q < 00, 

(37) tr (S|A)r /2 7„|A)| 1/2 ) - tr (5| A)| 1/2 7oo| A)| 1/2 ) , 

where l/p + l/q = I. 

We denote ||7„||f, p := || |-Do| 1 ^ 2 7n|-Do| 1 ^ 2 ||p- Given pi > P2 > 1, it is always pos- 
sible to extract a subsequence of 7n - again denoted by 7n -, such that it converges 
weakly in (F, ||.||f, P i) and (F, ||.||f, P2 )- Denote the weak limits respectively by 75? 
and 7^ 2) . Then for all B G & qi (£) we have B G 6 92 (ft) and 

tr(B\D \ 1 / 2 7 t ) \D \ 1/2 ) = lim tr( J B|Do| 1 / 2 7 n |^o| 1 / 2 )=tr(fl|D | 1/2 7^ ) |£'o| 1/2 ) . 

71 — *-00 

Thus, for all B G & qi (fi), we have 

tr( J B|^ | 1 /2 (7 (p 1 )_ 7 (p 2 ))| /Jo |i/2 ) = . 

This yields 7 &° =7^. 

Therefore, in the other parts of this proof, we will always assume we have chosen 
a subsequence of 7n such that, for the considered p, the weak-limits in (F, \\.\\f,p) 
exist and coincide. 

Step 2. Lower Semi-Continuity: We will now prove that taking the limit de- 
creases the energy, i.e. liminf„^ 00 f m (7„) > £ m (7oo). We define A m := X[ ,m](-D (,5) ) 
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and A' m := A+ — A m and split the energy functional. We treat the various terms 
separately 

£ m (jn) =Ti(j n ) + T 2 ( 7 „) + T 3 ( 7n ) + T 4 ( 7 „) + T 5 ( 7 „) 

:=tr(A^(£)W-m)A^ 7n )+tr {A m (D^ - m)A mln ) - ati (<^> 7lI ) 
+ atr(X^ 7n ) + aQ( 7 „, 7 „) . 
Step 2.1. Fix a basis (ee)e&i of each element being in H 1 (G). Then 
Ti(ln) := tv[A' m (D^-m)A' mln ) 

= tr ((A' m (D« - m)A; i ) 1 /2 7n( A; n (^(^) - m)A' m ) 1 /2) 
= 5> kl (A; n (^W - m)A^) 1 /2 7 „ (A ^(i)W - m)A^) 1 /2 efc ) 
* ' ' 

= £te (||X>o|- 1/2 /fc>(|i?o|- 1/a MPo| 1/2 7n|i?o| 1/a ) • 

k • ' 

-Hfc != 

Obviously, .ff/. is a non-negative Hilbert-Schmidt operator. Thus, applying first 
Fatou's lemma and then using (jSJJl, we get 

(38) 

liminf^ootr (A' m (D^ - m)A' m7 „) = liminf™ £ fe tr (fl" fc |A)| 1/a 7n|A>| 1/a ) 

> Efcliminfn^ootr^felDol 1 ^!^!!^) 

= E fc tr(F fe |D | 1 /2 7oo | Do |i/2) 
= tr {A' m (D^ ~m)A' mloo ) , 

which proves liminfn^oo Ti( 7 „) > Ti( 7oo ). 

Step 2.2. Because T 2 is continuous in the || • ||^ j2 -norm fLemma IA.13j) the 
claim follows for T 2 . 

Step 2.3. Since <p^ £ L 4 (R 3 ) (L emma lXjj) . we have, by [201 Theorem 4.1], 
that |Do| _1/ V (<5) l-Dor 1/2 € 64(8). By Holder inequality, this implies that T 3 is 
continuous in the || • |j^ 4/3 -norm, and lim n _> +00 T 3 (j n ) = T 3 ( 7 ). 
Step 2.4. We would like to prove 

(39) liminfQ[ 7rl ,7„] > Q[ 7 

oo ) Too J ■ 

n— *oo 

To that end, we will first show 

(40) lim j n (x,y) = 7oo(^,y) , 

n — >oo 

for a.e. (a;, y) £ G 2 , and 

(41) lim j n (x,x) = 7oo (ac, x) , 

n — >oo 

for a.e. x £ G. 

Now, ( 7n ) is a bounded sequence in 6 2 (,f)). Again, we can extract a subsequence 
such that 7n converges weakly to 7oo in & 2 {Sj). Using (|S7|l . we get for all B £ 6 2 (io) 

tr(B 7oo ) - lim tr(B 7n ) = lim tr ( \D Q \-$ B\D \~$ \D \^ n \D \^) 

ee 2 (f>) 

= tr (|D r*fl|A>r*|£o|*7oo|A>|') =tr(B 7oo ). 
Thus 7oo = 7oo . In particular we have 

(42) 7»(v)-»7oo(v), 
weakly in L 2 (G x G). 
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Using the spectral decomposition of the 7 ra 's, we may write each j n (x,y) as 

i 



ln{x,y) = ^uf\x)u) n) (y) , 



i=l 



where each sequence (u- ) nG N (i — 1, . . . , g) is an orthonormal family in £ i/ 1 / 2 (G)n 
A+^Sj. Since the sequence (tr (|p|7n))neN is bounded, it follows that, for each 
i £ {1, ...,q}, the sequence (t*j )neN is bounded in i? 1 / 2 (G). 

Therefore, applying O Theorem 16.1], for all \ £ C^(]R 3 ) and i S {1, . . . , q}, 

there exists a subsequence of (ttj )neN ~ also denoted (tq ) n eN ~ such that {x u ^) 
converges strongly in L 2 (G). 

Thus, after extraction of a subsequence of ( 

o I i denoted 

again by ( i4 , • • • j ) , we obtain, for all i = 1, . . . , q and for almost 

V / riGN 

every ieG, 

(43) ^(aO-^Or). 
Consequently, we obtain 

9 



(44) 7 n (x,y) -/3(x,y) := £ uj^fx) ^(y) , 



L 

i=l 

almost everywhere in G 2 . Now from (|42|l and (|44l) it follows by standard arguments 
that "foo{x, y) = [3{x,y) almost everywhere in GxG. Thus J n (x, y) converges almost 
everywhere to 7 0o (x,y). This together with 14411 implies (|40|) . 
The above also immediately implies that 

(45) 7oc=EK M )K-° 0) I' 

i=l 

proving ijlTjl. 

Applying Fatou's lemma to the pointwise positive functions 
(x,y) i — ► 7 n (x, x)7„(y, y) — |7„(x,y)| 2 and using in addition (|40|l and (|41|) yields 

liminfQ[7„,7n] > /liminf 2^3^ ~ ^i^ dy 

2 



f 7oo (x, x hoc (y, y - 7oo (x, y) 2 r , 

/ i i dxdy = Q 7oo,7oo 

J x-y 



which proves (|39ll . 

Step 2.5. Since t^ 5 ) £ L 4 (R 3 ) (Xemma [Q)) and < <^ 5 ) (by Lemma 
and the positivity of 1^ and 6), we have |£>W |-V2_x-(*) \D^ \-^ 2 £ 6 4 . Thus T 5 
is continuous in the || • || ^4/3-norm (by Holder inequality). Therefore 

lim T 5 (7„) = T 5 (7oo) ■ 



n — -+oc 



This concludes the proof of Theorem 13. 81 □ 
Theorem 3.9. Assume < S £ F, q £ N, q < Z and 

7ra(l/4 + max{tr S, q}) <(d — 4a tr 5). 
Then the functional £\„(s) (see the Definition (1231) and \ 1911 ) has a minimizer. 

Proof. Using lemma lA . 61 and the assumptions, we obtain 

(46) na< \ + max { tr l}) < c a,Z,m,S- 
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Therefore, by Theorem 13. 61 the functional £ m is bounded below on S q . 

Lemmata 13.11 13.31 and 13.41 together with Remark 13.51 imply that it suffices to 
minimize over positive 7's in Sq S ^ that can be written as a sum of a projection and 
a rank one operator. Moreover, Inequality (|46|l and Lemma |A. 121 permit to apply 
Lemma 13.71 which, together with the above and the fact that q is an integer yields 

(47) inf{£ m ( 7 )|7 G S<f>} = inf{£ m ( 7 )| 7 e S<J\ 7 2 = 7 }. 

Thus, applying Theorem 13.81 and again Lemma 13 . 71 shows that £ m \qW has a mini- 
ng 

mizer in So„ . Therefore the functional f Ls) also has a minimizer. □ 

4. Properties of the Minimizers: No-Pair Dirac-Fock Equations 

4.1. Infinitesimal Perturbations of Projections. In this section, we will prove 
that all minimizers 7 fulfill a no-pair Dirac-Fock equation. We first need to state 
preliminary results. The first was already used in the adiabatic theory (see Nenciu 
[T8"| and references therein). 

Lemma 4.1 (Nenciu). Given an orthogonal projection Pq, its orthogonal comple- 
ment P(J- = 1 — Pq, a bounded operator A, and £€l, with 4|e| \\A\\ < 1, there exists 
an operator B e with \\B e \\ < 4 ||A|| 2 such that 

(48) P e = P + e (P AP+- + P a ± A*P Q ) + e 2 B e 
is an orthogonal projection. 

Proof. We set 

(49) P e := / ± dz, 

2tti J| 2 -i|=i z - Pq - ea 

where 

a := P Q AP^ + P Q ± A*P . 

We observe that ||a|| < ||j4||. Therefore under the assumption 4e||vl|| < 1, we obtain 
that o-(P + ea) C (-1/4, 1/4) U (3/4, 5/4). Thus, by the holomorphic functional 
calculus, P c is the projector onto the eigenspace of (Pq + ea) corresponding to 
(3/4,5/4). 

e f 1 1 
P e ^p Q + —(j) — a — dz 

(50) 



27" J\z-l\=± Z ~ P Pq 



6 I 1 1 1 A 

® —a — a dz . 

27rl J\z-1\=% z ~ "0 z-Po-ea z-Pq 

Since Pq is an orthogonal projection, there exists a basis (ej)j^ of Sj and I C N 
such that 

Po=Tje j )(e j \. 



Note also for later purpose that 

(51) ^ = ^EN(^! E 1^-1 = ^0 + ^. 

jei jen\i 

Using Cauchy's residue Theorem and l|51|l. we obtain 

If 1 1,1/ /i in dz 
<P a dz = <f> (P(taP + P aPn) — r 

(52) 2znf lz _ ll= ^z-P z-Pq 2m jf[ 2 _ 1|= i V 1 z(z - 1) 

= P^taPo + PoaPn^ = a . 
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This proves that the second summand in the right hand side of (|50|) is equal to ea. 
This leads us to introduce : 

(53) B e :=—^<L — — a a — l —dz. 

Zm ,/|«— 1|=-| z — rj] z — Pq — ea z — Fq 

Since ||a|| < ||A|| and cr{P + ea) C (-1/4, 1/4) U (3/4, 5/4), we have 

1 1 1 A 

■ dz 



| z _l|=i z - P Q z- Po-ea z - P Q 

(54) < Trllall 2 sup ||^-|| 2 sup || i ||<16vrP|| 2 . 

| 2 _l|=i z- p o | z -i|=i z-P -ea 

This, together with l|50|l and (|52|l . proves Lemma |4~T1 □ 

In the case when 7 is an orthogonal projection with range in Sj+, we apply 
Lemma ETT1 with P = 7, P e — j e (given by iHHJl). 

Lemma 4.2. Assume < 5 F and 7 an orthogonal projection in Sg^ , q G N. 
Then, for operators A such that A^AA^ — A, \Dq\A 6 &i{Sj) and e sufficiently 
close to zero, 7 £ is again an orthogonal projection in S^j . 

Proof. By construction, since A^ AA^P = A, we have 

A^TeA^ - Af 7e A<? = A^ lt Af = . 

Moreover, 7 2 = 7e , thus -A^ < 7 e < and A^A^ = 0. 

The trace condition is obviously fulfilled since tr 7 e depends continuously on 
the parameter e. That £>o7e is trace class follows immediately from the explicit 
expressions for the difference j e — 7 in (|50|l . from 1|51|) and the assumptions on 
A. □ 

4.2. Minimizers are Projections. 

Theorem 4.3. Assume that < i5 £ F, g € N, and 

< 7ra(l/ 4 + max{tr <5, q}) < c a< z,m,b ■ 
If '7 is a minimizer of £L(s>, £/ien 7 = 7* = 7 2 = A+ 7AY 5 . 

Proof. The proof of 7A+ = 7 is a consequence of Lemma f3.ll The proof of 
7 2 = 7 follows exactly the lines of Lemma 13.41 except that the iteration of the 
process is superfluous here. □ 

4.3. Minimizers Fulfill the No-Pair Dirac-Fock Equations. Eventually we 
derive the Euler equations for the minimizer of the energy. 

Theorem 4.4. Assume q £ No and 7 is an orthogonal projection minimizing £ in 
S'dq ■ Then 7 commutes with the no-pair Dirac-Fock operator A+D^A^, i.e., 

(55) [ 7 ,Af D^Af ] =0. 

Proof. Let A S B(Sj) such that 

(56) D^A* e B(Sj) . 
Then, for e sufficiently close to zero, the projector 

(57) 7e := 7 + ea + , 
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with 

(58) a = 7 A$?i4AjV + 7 ± Af A*kf 7 , 

and B c given by l(53"|) (with P Q replaced by 7), belongs to Sg q '. Moreover 

(59) £(7 e )-£(>y) 

= e{tr (D z a) + 2a^Q[j, a]} + e 2 {tr {D z B e ) + 2a s RQ["f, B e ]+aQ[a+eB e , a+eB e ]}. 

We want to show that the last term is o(e). By 1)13(1 . 114(1 and Lemmata IA. 71 and 
IA.8I it is sufficient to show that there exists a constant c < 00 such that, for all 

ee (-1,1), 

(60) max{||B £ ||i, ||a||i, P^lli, P (5) a||i} < c . 
We first have 

^ IpWaUr < ||DW 7 Af ^Vlli + P ( VA^A*Af 7II1 
(61) 

< ||^ 7 ||i \\A\\ + \\D^\\ \\A*\\ H7II1 + \\D&A*\\ || 7 ||i < c , 
since j £ F and ((56(1 is assumed. We also have 

|p (5) iyi = i|| / DW— a a— dz||! 

2tt /| z _i|=i 2-7 z-7-ea z - 7 

* f />> IT tV-I- »T-t^ "TTT^H d " 

27rJ l + ie^-7 l + ~e^- 7 -ea 1 + ±e l f - 7 

<ci||o|| / ||L> (l5) f A 7 + 1 1 — (1 -7)1 a||i<ty> 

" " " 7 Vl + ie^-l' l+|e^ v "J 111 r 

< 02 ||a|| (||# (5) 7l|iN| + ||# (5) a||i) , 

where the constant C2 is uniform in e for e close to zero. We also have used (15 111 . 
Using ((61(1 . and 7 G F yields 

(62) \\D^B e \\<c. 
Similarly to the above, we show 

(63) N|i<c and ||S e ||i<c. 

Inequalities 03} , JUS and (JHSJ yield tfHHft . 

Since £( 7c ) — £(7) > 0, whatever the sign of e is, we conclude that the term 
linear in e in ((59(1 has to vanish 

(64) tr (D z a) + 2a5RQ[ 7 , a] = tr (D^a) = . 

Thus, for all A satisfying equalities JHEJ) and JHU) and the fact that [Ay, 7] = 
(since 7 is an ortho normal projection in Sgj) yield 

tr (^A^D^AfjA) + tr {jA^D^A^^A*) = . 

Replacing A by iA, we first obtain : 

tr( 7 x Af £>«a£> 7 A) = tr ( 7 A^ D^A^ 7 X A*) = . 

Since A can be taken in the set of rank one operators of the form \u) (v\, with u 
and v in C^°(IR 3 ) <g> C 4 , we obtain 

7 x Af flWA V7 - -yAyWAyV = , 

which yields l(55|) . □ 

This result can be also written as follows 
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Theorem 4.5. Assume that 7 is an orthogonal projection minimizing £ in Sg q . 
Then there exist q normalized spinors ipi, G A+\S}) n ®(D Z ) such that 

(65) 7 = i><XlM, 

i=l 

and 

(66) AjWAf^ = e,-^, for i = 1, . . . , g , 
wifft ei, e q G [0, 1]. 

Proof. The proof is immediate since the range of 7 is finite dimensional reducing 
it to the simultaneous diagonalization of two commuting Hcrmitian matrices. □ 

One may characterize the eigenvalues ei , . . . , e q more precisely following Bach et 
al 0: 

Theorem 4.6. [There are no unfilled shells in no-pair Dirac-Fock theory] Un- 
der the same assumptions of Theorem \4-5\ t\,...,e q are the q lowest eigenvalues 
of A^ . Moreover, if e q +\ denotes the (q + l)-th eigenvalue (counting 

multiplicities) of the no-pair Dirac-Fock operator AyDWAr, we have, for all 
i = 1, •• ■ , q, the strict inequality ej < e q +\. 

Proof. Assume by contradiction that there exists a normalized eigenspinor v of 
with eigenvalue e such that e < maxjei, ...,e 9 } and not in the range 

of 7. Then, for a normalized eigenvector u of 7 with (u, A+D^A+u) > e and for 
7' := 7 — \u) (u\ + \v) (v\, we have 

£(*)(7')-£W( 7 ) 

= (v, D z v) - (u, L» z u> + 2o*Q [7, |«) («| - |u) («|] 

+Q[|«) («| - |«) H, |«) («| - |«) («|] 
< («,I)Wu) - (u,D { ' f} u) < 0, 

where - as in the proof of Lemma [3.41 - we used l|36l) . with S — \v) (v\ — \u) (u\. 
This gives a contradiction to the property that 7 is the minimizer. □ 

5. Outlook 

In this final section, we first express the energy as a functional of the unrenor- 
malized density matrix T and the Dirac sea A_. This has the advantage that the 
dependence of the energy on A_ becomes explicit and that the constraining condi- 
tion — A_ < 7 < A + becomes < T < 1, i.e. independent of A_. Throughout this 
section we will use the notation A := A_. 

For a given q E N, let us define the set 

T q := j(r, A) e 2$(f)) 2 I r, A orth. proj., (r — A) G Si, tr (r — A) = q, 
D Z (T - A) e 61, [D Z ,A] G B{Sj), [D Z ,T] G B(i5)} 
and the following functional on T q : 

g(T, A) := tr ((D z - m)(T - A)) + aQ[r — A, T — A] , 

where Q{-, •] is defined in l|24[l. Note that if A = A' 5 ' for some S G F, and if 
(r - A) G S 1 ^'), then wc have g(T, A) = £(T - A). 
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Theorem 5.1. Assume that (r, A) £ T q is a critical point of (£. Then, with 

7-r-A, 

(67) [D™,T\ = [D™,A] =0. 

Proof. For all e € K and A € such that : 

(68) 

D Z A G 6i(i5) , £>zA* G 6i(fl) , 4|e|||A|| < 1 , 8ae 2 ||A|| 2 < 1 and e 2 p||i < 1 , 
we define (see Lemma f4. II) the orthogonal projector : 

(69) A e = — !— <£ i dz , 

2m J| 2 _ 1 | = i z — A — ea 

where a = AAA 1 - + A 1 - A* A. From Lemma f4. II we get the decomposition 

(70) A £ = A + ea + e 2 B e , 

with Po replaced by A. 
Let us prove that (r, A e ) belongs to T q . 
We first show V - A e e &i(Sj). We have 

(71) ||o||i < 2||A||i < oo , 
and, as in the proof of (|54(l . we get 

(72) H^Hi < 167r|m| IHIi . 

Therefore, r-A £ = r- A + ea + e 2 B £ G 

We next establish D z A e G 6i(i)). Since (r,A) G T q , [D Z ,A] is a bounded 
operator and A, D Z A G &i(Sj). Thus 

H^AAA^Hi < WADzAA^U + \\[D Z , A]AA ± \\ 1 
<\\D Z A\\ 1 + \\[D Z ,A]\\ ||A||i<oo. 

Similarly, we can prove 

(74) ||D Z A X A*A|| 1 <||DzAli + ||[D^, A x ]|| ||A*||i<oo, 
which implies, together with i|73|) that there exists a constant c such that 

(75) \\D z a\\i<c. 

Using again that [ D z , A ] is bounded and Formula (|51Jl , valid with A instead of Pq 
and A ± instead of Pq , we get the existence of a constant c such that for all e small 
enough 

(76) WDzB^Kc. 

Inequalities (JZHJ and (JZEJ) yield D z A t G &i{Sj). 

Now, from (|7K|). i fTHjl and [Dz,A e ] = [L>z, A] + [£> z , ea + e 2 B e ), we obtain 
[D z , A e ]eB(Sj). 

We finally prove that tr (r — A e ) = q. For that purpose, we first note that, due 
to Effros ( see a l so Avron et al U Theorem 4.1]), and since from the above 
r — A £ G 61(f)), and both T and A £ are projections, we have tr (r — A £ ) G Z. 
Furthermore, from (|Tf)|) . 171|l and Q72J1 we get 

||(r-A)-(r-A e )|| 1 = 0(e). 

Since tr (r — A £ ) is an integer and tr (r — A) = q, this yields, for e small enough, 
tr (r - A £ ) = q. This concludes the proof that (r, A e ) G T q . 
Let us now prove that for 7 := T — A we have 

(77) [D M ,A] = . 
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Since (T, A) is a critical point of £, for all A 6 &i(Sj) satisfying l|68[) . we have 

<9£(r,A £ ) 



= 

e=0 



< 78 > Qe 
where A e is defined by (ZD} - ^ n ^ ne other hand, we have 
(79) 

£(I\A e ) =tr ((^z-m) 7 ) + aQ[7, 7 ]-e{tr ((£> z - m)o) + 2a3?Q[a, 7 ] } 
+ e 2 {tr {{D z - m)B e ) + aQ[a + eB t ,a + e£ £ ] - 2a5RQ[ 7 , 5 e ]} . 

Inequalities l|72|) and (|76|1 imply that there exists a constant c\ such that for all e 
small enough 

(80) |tr {{D z - m)B e )\ < c x . 

Furthermore, using Lemma I A. II we have 

tr(|A>| 1 / 2 |a + ei? e ||A ) | 1 / 2 ) = trda + e^l^l^oHa + e^l 1 / 2 ) 

< -tv^a + eB^DzWa + eB^ 2 ) 
a 

< i|||I> z ||o + eB 6 |||i - -\\D z {a + eB,)\\x, 
a a 

which yields, together with l|75|l and 176JI . a + eB e S F. Thus Lemma \'2 . 21 implies 
(81) 

\Q[a + eB e ,a + eB e }\ < 2D[a + eB e , a + eB e ] < |||a + eB e \\ i tr {\D \\a + eB e \) 

< ^Ha + eSellitr {\D z \\a+ eB e \) 

< ^||a + eB e ||i||D z (a + eB e )Hi 

According to i|75|l and i|76|l . we conclude from 1|81|) that there exists a constant C2 
such that for all e small enough 

(82) \Q[a + eB tl a + eB e ]\ <c 2 . 

Now we prove that there exists C3 such that for all e small enough 

(83) \Qky,B e ]\<ca. 

We have, for W^ 3 ^ being the mean field potential associated with J5 e S F, as 
defined in Q| 

\Q[-y,B e ]\ = \tr{W^ 7 )\- 

Moreover 

(84) |tr(W^> 7 )| < ||W^7l|i < W^iD^W WD^Dz^W \\Dzl\\i • 
Using Lemmata I A . 31 and I A . 41 with B e instead of <5 implies 

(85) WW^'HDo)-^ < 4115.111 • 

Since Dq{Dz)~ 1 is bounded, by using Inequalities Ij84(l . (|85|) and the fact that 
(r, A) g T 9 , we obtain lfS5)). 

Collecting JE3, and (jHSJ yields, together with f75^. 

£(T, A e ) = tr {{D z - m) 7 ) +aQ[ 7 , 7] -e{ tr ((D^ - m)o) + 2a5JQ[a, 7] } + 0{e 2 ) 

= g(T, A) — etr ((£>W - m)o) + 0(e 2 ) 

= g(T, A) - e tr (A x (£>^ - m)AA) - e tr {A{D^ - m)A ± A*) + C(e 2 ) . 
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For A fixed as above, this result remains true for all e small enough. Therefore, 
(JZHJ implies 

tr (V(L>^ - m)AA) + tr (A(D™ - m)A ± A*) = . 

As at the end of the proof of Theorem 14.41 we obtain [D^\ A] = 0. 

Finally, exchanging the roles of T and A in the above proof yields [D^\ T]. □ 

Appendix A. Some Spectral Properties of Screened Coulomb-Dirac 

Operators 

We recall the following result: 
Lemma A.l (Brummelhuis et al 8 ). Let 

d := (1/3) (1 - (aZ) 2 ) 1/2 ((4(a^) 2 + 9) 1/2 - AaZ) 
and assume 0<aZ< y/3/2. Then 

(86) \D Z \ 2 >d 2 \D \ 2 and \D Z \ > d\D \. 
We would like to compare \D \ and \D^\. 

Definition A. 2. Given positive a, Z , and m and 5 E F , we define 

(87) c Q)W := sup{c S R | \D^\ > c\D \} . 
Lemma A. 3. If 6 E F, then, for all u E ii^R 3 ) ® C 4 , we have 

\\X {S) u\\ <2||5||i||Vu|| . 
Proof. Using the spectral decomposition of <5, we have 

\\X^uf=l 6{X \ )W f )u{z) ^d X dydz 
J |x-z||x-y| 

fi J l x - z ll x -y| 

<£/ WIWiMII&MI 

<4^ J lAiMAjl Mx)\ Mx)\\\VufAx 
<M\6\\l\\Vu\\ 2 , 

where we used the Cauchy-Schwarz inequality, the Hardy inequality, and the iden- 

tityEi MHHIi- □ 

Similar estimates can be found for the direct part, 
Lemma A.4. If S E F, then \\^ s) u\\ < 2 ||5||i||Vu||, for all u E H 1 {M. 3 )®€ i . 



|x-y| 



■Ay 



dx 
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Proof. We have successively 



*(■>■■■•■)»{:. :.) ( / , H ff , dy)<Uh 
x-y |z-y 



< A J \5{x,x)\Ax J |d(*,*)|ds||V«|r<4||«||i||V«|r, 

where we used Hardy's inequality in JSSJ. □ 

A direct consequence of Lemma IA.4I and the fact that square root is operator 
monotone is 

Lemma A.5. If SeF, then \cp^\ < 2||<5||i|A,| • 

Lemma A. 6. If S G F and aZ < y3/2, the following operator inequality holds 

(89) > (d-4a\\8\\i) . 

Proof. This is a direct consequence of Lemmata lA.3l and IA.4I since we have for all 
MinSJpW) 

\\D^u\\ = || (p z + a^ 5) - al (5) ) u|| 

( 90 ) > \\D z u\\ - a\\tp {5) u\\ - a\\X {5) u\\ 

> (d-4a||J||i)||I>o«|| . 

Therefore |DW| 2 > id — 4a||<5||i) 2 |Z)o| 2 , which concludes the proof since the square 
root is operator monotone. □ 

Lemma A. 7. Assume SeF ande>0, then pW, <p(l ff D e L 3+e (R 3 ) n L°°(R 3 ). 

Proof. As before we denote by (A„)„ g n the eigenvalues of S and by (£„) a cor- 
responding orthonormal basis of eigenfunctions. Since |<^^^|(x)| < tpW)(x), it is 
sufficient to prove the result for ip^ s ^ . 

We first prove that ip" s " e i 3+e (R 3 ). We write XR f° r the characteristic function 
of the ball with center and radius R, and set Xr '■= 1 — XR- We get 

\\\Cn\ 2 * | • |- 1 || 3+e < || |y 2 ||f || XR I • ll3/(2 + e) + || XR I • I" 1 lU+e 

<c[<e„, |pie») + i] , 

where we used the Hausdorff- Young inequality and the Sobolev inequality. Multi- 
plication by |A n |, summation over n, and the triangular inequality yields the result. 

\W m h + , < , |V| U + 1] < C\\5\\ F . 

n 

Next we estimate H^'^Hoo- Using Kato's inequality, we get 



»(x)| <^|A n | J dy\Mv)\V\x-y\ < \ El A «l^' IpI&> 



□ 

Lemma A.8. If SeF, then \X^ \ < ip^ . 

Proof. This is a straightforward consequence of the spectral decomposition of S and 
the Cauchy-Schwarz inequality. □ 

Lemma A. 9. IfSeF, then is relatively compact with respect to Dq. 
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Proof. Since tp^ € L 4 (M 3 ) by Lemma TA. 71 using an inequality of Seiler and Simon 
fM Theorem 4.1], we have 

(91) M 5) (-A + m 2 )- 1 / 2 ^ < ||^ (5) (-)ll4 ||(| • I 2 + m 2 r 1/2 |U < oo , 

implying £>o~ V W e 6 W . 

We next prove that X^ is relatively compact with respect to Dq. Let us denote 
by <5_|_ (respectively <5_) the positive part (respectively the negative part) of 5: 
5 = 5+- 5-, 5± > 0. 



We have, using 5±(x,y) = 5±(y,x), 



5±(x,y)8±(y,z) 



(u,(X^) 2 u) = I u(x) (J "^°'"=™'~{ u(z)dydz ) dx 
I ^rG{x,z)iG(z,x)i^-dxdz, 



where G[x, z) 



u{x)\ and h is positive and measurable. 



x-yi |y- z l 

Applying the Cauchy-Schwarz inequality and |<5±(x, y)\ 2 < S±(x,x)6±(y,y) yield 

h(x) 2 f 5±{x 1 y)5±{y,z) 



{u,(X {S±) ) 2 u) < 



h(z) 




-dy 



|u(a;)| 2 da;dz 



Picking h{x) = 5±{x,x) 1//4 implies 
(92) 

Mx«*>)'u> < / Hxf (/ l^fj^ v) 



\u(x\ 2 (J ^fe)(y)MMldj/) dx < M^IU^M 



Now, similarly to l(9*T|). we have ||(/3^ 5± '(— A + m 2 ) 1//2 ||4 < oo. Thus, since -D 1 
is bounded, we obtain Dq 1 <p {S±) Dq 1 e & 4 (Sj). Mor eover, by Lemma fA. 71 we have 
||<^ (5±) ||oo < oo. Therefore, using |g2J, we get D 1 (X^) 2 Dq 1 E 6 4 (ij) and thus 
X^Dq 1 e & 8 (Sj) C Soo(J3), which implies X^Dq 1 compact. Using this, and 
the fact that is relatively Do compact conclude the proof. □ 

Lemma [A. 91 the Kato-Rellich Theorem, and Weyl's theorem imply 

Lemma A. 10. Let Z > such that aZ < 1 and let 5 € F. Then the operator 
£)W = Dz+a^'-X^) is self-adjoint with domam r D(D^) = D(D Z ). Moreover 

(Jess (-D (<5) ) = (T oss (D Z ) = ( CO; -m] U [to, +oo) . 

If (5 is a positive density matrix of q electrons with q < Z + 1, then the potential 
of the nucleus prevails giving an attractive Coulomb tail at infinity. This leads us 
to expect that has infinitely many bound states in the gap accumulating at 
m. The following theorem expresses this expectation formally. 

Theorem A. 11. Assume < 5 S F and tr 5 < Z. Then the operator £)W /ias 
infinitely many eigenvalues in (0, to) accumulating at the point to. 

Proof. Pick a function / € Cg^R 3 ), normalized in L 2 (R 3 ), which is also spherically 
symmetric. We also define, for all R > 1, the functions 

f R (x) := R- 3 / 2 f(x/R) and ^ = (/«, 0,0,0) . 
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Straightforward calculations using scaling arguments and the specific form of i/jr 
give 

(93) 

|| (D -^ + a^-aX^i; R \\ 2 

= \\a- ivV^f + \\m[3M\ 2 + II (-^ + ^ {S) - M\ 2 
+ 25? (mPi/jR, a ■ tW>h\ + 23? / (~r^ + aL P (5) - ip R) a ■ ^Vi/jr 

+ 25? 



(m^R, (~ + - crfOj 
= || V^|| 2 + m 2 + a 2 || (-|r + <p& I< 5 ») <fo|| 

+ 2ma /fe (-]|f + ^ - Xi5) ) 
< ^l|V/|| 2 + m 2 + ^(||^/|| 2 + 8( ? 2 ||V/|| 2 



+ 2m 



i c i / ( 

< m + — + 2m (ipR ' 



x 



where ci is a constant independent of R. Note that we have used Lemmata IA.3I 
and I A. 41 for getting an upper bound of || (^ - X^)ip R \\ 2 . 

Now, since / is spherically symmetric, using Riesz's rearrangement inequality 
and Newton's Theorem, we obtain 



^Hn) = / lf f )l P f ] dxdy < / MM dxdy 
/ 7 |x-y| V |x-y| 

/• |/ fl (x)| 2 /■ /• |/ fl (x)| 2 

< J ^^tej Ps (y)dy< q J ^dx, 



where we used, since p^(y) = X)t=i 2/) * s P°siti ve j that 

pKy)dy=||p|||i = ||p 5 || 1 = g / . 



Note also that (ip R ,X {s ^ R ) > and 



with C2 := (/, t^t/) > 0. Thus, we have 

c 2 (Z - q') 



(94) 



R 
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Now dnSJ and (§D imply 

|| [D&-E] Vh|| 2 =||-D^fl|| 2 + J B 2 -2^(^,i)WVR) 

< m 2 + ^ + 2m<^„ a (-ft + ^ - X^) Vii) + £ 2 

-2£<^, a (-ft + tpW - !(*)) Vfl) - 2S(^, A)V>fi) 

< m 2 + + 2(m - a (-^ + - A< 5 )) ^) 
+E 2 + §r - 2mE 

= (m - £) 2 + 2fr - 2(m - E)a C2(z R q " > . 
Therefore, if R is large enough, we get the inequality 

|| (d^-e) Vflll <\m-E\, 

which implies, by taking E — m/2, that D^ s > has at least one eigenvalue Ai in 
(0,m). 

Now by taking E = (m + Ai)/2, by the same argument as above, one gets a 
second eigenvalue A 2 G (Ai,m). The iteration of this procedure yields an infinite 
sequence of eigenvalues (A„) of D^ s ' in (0,m) tending to m. □ 

A similar result holds for the no-pair Dirac-Fock operator. 

Lemma A. 12. Assume that 5 and 7 are two positive definite finite rank density 
matrices. Assume, in addition, that 7 is purely electronic having particle number 
not exceeding Z, i.e., we have 

SeF,5>0,-/eFr\S s , A^A^ = 7 and tr 7 < Z . 
Moreover we assume 

d- 2a(2tr 5 + tr 7) > . 

Then the no-pair Dirac-Fock Hamiltonian A^P A^ has infinitely many eigen- 
values in (0, m). 

Proof. We first prove that 

(95) <WAj } £>«Af ) = [m, +00) 

Lemma [A . 91 implies that VF^-D^ 1 and W^ 7 - 1 /?^ 1 are compact. By assumption, we 
have (d- 4a||<5||i) > 0; thus JU) implies that DqCD^) -1 is bounded. This yields 

and 

A f ( W « - ^W)A^(Af dWaS?)" 1 g Sooft) . 
Using the Kato-Rellich Theorem, this implies 

a css (A^ 7 >A^) = ^(Af (DO + W« - W^)Af) = a css (^D^A^) . 

Together with Lemma [A. 101 this proves 

Set q 1 :— tr <5 and q :— tr 7. We denote by tp^ and X^' respectively the direct 
and exchange operators associated to 7, defined by replacing S with 7 in 10 and 
©. For all u G Af? (H^R 3 ) <g> C 4 ), we have 

(u, L> (7) u) = (u, £> (,5) u) + a(u, <p M u) - a(u, X {l) u) - a(u, (p {5) u) + a(u, X {5) u) 

(96) < -{u, |D (5) |u> + a(w^ (7) «> > 

where, in Inequality IjHfiJl . we used Lemma fA.8l and the fact that A^ 7 ' > 0. Now, 
from Lemmata IA.5I and IA.6l we obtain 

(u, D (7) u) < (-(d - Aaq') + 2aq) (u, \D \u) . 
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Since — (d — Aaq') + 2aq < 0, we have 

(97) A^WA? < . 

Now, thanks to Ij97(l . one can apply 0| Theorem 3]. With their notations, since 
D(D^) = CpW), we first define 0± := f^(D^j\ n A^ij. Then, denoting by 
p n (A) the n th lowest eigenvalue of the operator A, we obtain 

m>n n (DW\. Msi ) > inf sup ((, D^() . 

v 7 M + ca + ceM+ ea_ 

dim(A'/+)— n 

Moreover, 

inf sup (C, C (7) C) > inf sup ((, D^\) = /x„ (a^D^A^) • 

dim(M + )=n ||f|| = l dim(M+)=n ||£||=1 

Therefore 

m > M „ (s W | A w«) > Mn (Af DWAf ) . 

Since from Theorem IA.11I the operator has infinitely many eigenvalues in 

(0, to), using J^SJ, we finally get the expected result. □ 

If we have a non-negative spherical symmetric density p with 




p{y)dy < Z , 



then 

-z/h+p* i - r 1 > o . 

This implies that the n-th eigenvalue of + </? can be estimated from below by the 
n-th eigenvalue of Dz-q- For a non-spherical symmetric potential this situation for 
the positive eigenvalues is disturbed only slightly. 

Lemma A.13. Assume < 6 e F. Then X(Q, m ){D [S) ){D {5) - to) is a Hilbert- 
Schmidt operator. 

Proof. Since 6 € F, by writing W^D Z X = W^D^DoD' 1 , Lemmata FOl and 
IA.9I imply that is relatively compact with respect to Dz- Moreover, since 

5 > 0, we have > (see lemma FO|l . 

Let Ao(0) < Ai(0) < ... be the ordered positive eigenvalues of Dz, including 
multiplicity. We first prove that for all e £ [0, 1], there exist N+(e) and AT_(e) in 
{-oo} UZU{+oo}, with N-(e) < N+(e), and (Afe(e)) fee{Ar _ (e) jAr+(e)} such that 

i) a (D z + eaW^) n (-to, m) = {A Ar _ (e) (e), A Ar _ (e)+1 (e) ! . . . , Ajv + ( e )(e)} 

ii) AAr_( e) (e) < AAr_( £)+ i(e) < ... < A w+ ( £ )(e) 

iii) For all k € {./V_(e), . . . , iV + (e)} the functions Afc(e) are continuous, monotone 
increasing. 

Using Kato's perturbation theory for isolated eigenvalues and numbering the 
eigenvalues that are crossing with respect to their ordering - namely the largest after 
crossing gets the highest index -, yields i) and ii). Continuity is also a consequence 
of perturbation theory for eigenvalues. The asserted monotonicity is a consequence 
of positivity of . 

We next show that the number M of eigenvalues A&(e) that crosses zero when e 
increases from to 1 is finite. By the Birman-Schwinger principle, 

M = #{A e [1, oo) | A eigenvalue of R} , 
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where R = — (W^) 1 ^ 2 D Z 1 (W^) 1 ^ 2 is the Birman-Schwinger kernel with energy 
0. Thus, since we count only the eigenvalues of R larger than 1, we have, using the 
notation W := W {s) and ip := f (s \ 

M < \\R\\\ < WW^ip-^ip^lDol-ilDo^Dg^Dol^Do^tpiip-WWi 



(98) 

< \\W^-^\\ 8 || lA^ZVH 8 ll¥>*lA>| 
Since < = f (s) - < </jW, we get 

(99) ||(^ (5) ) 1/2 ((^ (5) )- 1/2 || < 1. 
Lemma IA.1I yields 

(100) |||A)| 1/2 ^ 1/2 || <d^ <oo. 

Using H86|) and (501 Theorem 4.1] as for Inequality 1|91|) in the proof of Lemma 
yields 

(101) iK^v^oi-v^ < ||^)(.)||i /2 ||VI-l 2 + - 2 llI /2 < oo. 

Collecting P fy -fiTE Q proves M < oo. 

Therefore, apart from the eigenvalues Ao(l), Ai (1), . . ., the operator has only 
finitely many other eigenvalues in [0, m]. Thus, X(o,m){^ ){^ — m) e &2(Sj) 
follows if the series Sk>o(Afc(l) — m ) 2 — S/c>o(^ fc (0) — TO ) 2 ^ s convergent. At this 
point we remind the reader that the relativistic hydrogen eigenvalues A& (0) can be 
grouped into "multiplets" of 2n 2 eigenvalues corresponding to one non-relativistic 
eigenvalue m — Z 2 a 2 /n 2 . Each element of such a multiplet can be bounded from 
below by the previous non-relativistic eigenvalue m — Z 2 a 2 /(n — l) 2 . Thus, up to 
an unessential multiplicative constant, X) n >2 n2 /( n ~ I) 4 is a convergent majorant. 
This proves the claim. □ 
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